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Abstract 

A detailed investigation of the energy loss of gluons that traverse a thermal gluonic medium 

simulated within the perturbative QCD-based transport model BAMPS (a Boltzmann approach 

P-f to multiparton scatterings) is presented in the first part of this work. For simplicity the medium 

P-T 

q_) ' response is neglected in these calculations. The energy loss from purely elastic interactions is 

4= : 

compared to the case where radiative processes are consistently included based on the matrix 
element by Gunion and Bertsch. From this comparison gluon multiplication processes gg — > ggg 
qq | are found to be the dominant source of energy loss within the approach employed here. 

■^J- ■ The consequences for the quenching of gluons with high transverse momentum in fully dynamic 
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simulations of Au+Au collisions at the RHIC energy of y/s = 200 AGeV are discussed in the second 
major part of this work. The results for central collisions as discussed in a previous publication are 
revisited and first results on the nuclear modification factor Raa for non-central Au + Au collisions 
are presented. They show a decreased quenching compared to central collisions while retaining the 
same shape. The investigation of the elliptic flow V2 is extended up to non-thermal transverse 
momenta of 10 GeV, exhibiting a maximum V2 at roughly 4 to 5 GeV and a subsequent decrease. 

Finally the sensitivity of the aforementioned results on the specific implementation of the effective 
modeling of the Landau-Pomeranchuk-Migdal (LPM) effect via a formation time based cut-off is 
explored. 

PACS numbers: 12.38.Mh,24.10.Lx,24.85.+p,25.75.-q 
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I. INTRODUCTION 



It has been established by experiments at the Relativistic Heavy Ion Collider (RHIC) that 
jets with high transverse momenta are suppressed in Au + Au collisions with respect to a 
scaled p+p reference l|,|2j. This phenomenon of jet quenching [3| is commonly attributed to 
energy loss on the partonic level as the jets produced in initial hard interactions traverse the 
hot medium, the quark-gluon plasma (QGP), created in the early stages of such extremely 
violent heavy ion collision. Due to the large momentum scales involved the energy loss 
of partonic jets can be treated on grounds of perturbative QCD (pQCD). Most theoretical 
schemes attribute the main contribution to partonic energy loss to radiative processes, where 
gluons are emitted in bremsstrahlungs-like interactions |4H10|. 

It is a major challenge to combine these approaches with models that describe the bulk 
properties of the medium, most notably the strong elliptic flow, quantified by the Fourier 
parameter v-i- Com par ison to hydrodynamical calculations shows that the viscosity of the 
QGP is quite small [11] . possibly close to the conjectured lower bound 2 = from a corre- 
spondence between conformal field theory and string theory in an Anti-de-Sitter space 12]. 
Results from hydrodynamical simulations are used as an input for the medium evolution 
in jet-quenching calculations (see [k| for an overview) and as ingredients in sophisticated 
Monte Carlo event generators 14J. However, these approaches treat medium physics and 
jet physics in the QGP on very different grounds. In a previous publication |15j] we have es- 
tablished that partonic transport models including radiative pQCD interactions can provide 
means to investigate bulk properties of the QGP and the evolution of high-energy gluon 
jets within a common physical framework. 

There the build up of elliptic flow and the quenching of gluon jets have been studied in 
fully dynamic microscopic simulations of a gluonic medium created in Au + Au collisions 
employing mini-jet initial conditions. Matching the experimental values on integrated V2, the 
nuclear modification factor Raa was found to be completely flat and slightly below the value 
obtained by Wicks et al. in their calculations for gluonic Raa [lOi] - We clearly demonstrated 
that the incorporation of inelastic bremsstrahlung processes into partonic transport models 
provides a promising approach and that the consequences of such a description need to be 
investigated carefully. 

In this work we focus on the energy loss mechanisms within our parton cascade BAMPS 
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(a Boltzmann approach to multiparton scatterings) [lfj, [l7j and in particular provide details 
on the evolution of high-energy gluons that traverse a thermal gluonic medium. We present 
our results in such ways that will allow other partonic transport models to easily compare 



with our findings |36j. Note that the term "jet" as used throughout this paper refers to a 
single parton with high energy and does not fully coincide with the experimental notion. 

This paper is structured as followed: In section [XT] we introduce the transport model 
BAMPS used in our investigations. We especially focus on the details concerning the inelastic 
gluon multiplication processes and the involved effective implementation of the Landau- 
Pomeranchuk-Migdal (LPM) effect. In section UTT1 we discuss the evolution and energy loss of 
high-energy gluons in a static medium exclusively considering binary gg -H- gg interactions. 
After having established this baseline, in section HVl we turn towards the evolution of gluonic 
jets in a static medium including inelastic gg <-> ggg processes. Afterward, in section [V] 
we extend our investigations of central Au+Au collisions at RHIC energies from a previous 
publication to non-central collisions with an impact parameter 6=7 fm. In section I VI I we 
investigate the sensitivity of the previously obtained results on details of the modeling of 
the LPM effect via a cut-off in momentum-space. 

II. THE TRANSPORT MODEL BAMPS 

BAMPS [if], [l?! is a microscopic transport model aimed at simulating the early stage 
of heavy ion collisions on the partonic level via pQCD interactions consistently including 
parton creation and annihilation processes. The use of microscopic transport calculations 
provides a realistic way of understanding the complex evolution of the QGP including the 
full dynamics of the system and allowing for non-thermal initial conditions. Partons within 
BAMPS are treated as semi-classical and massless Boltzmann particles. At this stage the 
model is limited to gluonic degrees of freedom though the implementation of light quarks is 
underway and will be presented in an upcoming publication. The investigations presented in 
this work, however, are exclusively dealing with high energy gluons traversing purely gluonic 
media and thus Nf = is understood for all calculations. 

The interactions between partons are implemented based on leading order pQCD matrix 
elements from which transition probabilities are computed. These are used to sample the 



interactions of particles in a stochastic manner 16j . The test particle method is introduced 
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to reduce statistical fluctuations and is implemented such that the mean free path is left 
invariant. For elastic interactions of gluons, gg -H- gg, we use the Debye screened cross 
section in small angle approximation 

dVgg^gg _ ^ 7ra s /i\ 

dql (qi + mlf ' 1 ) 

Throughout this work we employ a fixed coupling constant a s = 0.3. The Debye screening 
mass is computed dynamically from the local particle distribution / = f(p,x,t) via 

/d 3 p 1 
j2nyp Ncfj (2) 

where da = 16 is the gluon degeneracy factor for N c — 3. 

Inelastic gg -H- ggg processes are treated via an effective matrix element based on the 
work by Gunion and Bertsch 18]. Detailed balance between gluon multiplication and an- 
nihilation processes is ensured by the relation \M- g g^gg g \ 2 = do \M ggg ^ gg \ 2 . For the case of 
bremsstrahlung-like processes the matrix element employed in BAMPS reads 

lu ,2_ 72n 2 ay 487ra,qj _ 

\Mgg^ggg\ ~ (q 2 + ^2 k 2 [(k± _ q±)2 + ^ u ^9 V ■ 

qx and denote the perpendicular components of the momentum transfer and of the 
radiated gluon momentum in the center of momentum (CM) frame of the colliding particles, 
respectively. 

When considering bremsstrahlung processes the LPM-effect a coherence effect 

named after Landau, Pomeranchuk and Migdal, needs to be taken into account that leads to 
a suppression of the emission rate for high-energy particles. Since such an interference effect 
cannot be incorporated directly into a semi-classical microscopic transport model such as 
BAMPS, we choose an effective approach by introducing the Theta function in ([3]). This 
implies that the formation time r of the emitted gluon must not exceed the mean free path 
of the parent gluon A g , ensuring that successive radiative processes are independent of each 
other. See 20J for a recently developed probabilistic algorithm to include the LPM effect 
in Monte Carlo event generators. The transfer of this approach to a full-scale microscopic 
transport model, however, would pose an enormous numerical challenge. 

When comparing the formation time to the mean free path, i.e. the time between suc- 
cessive interactions of the parent gluon, special attention needs to be paid to the frames of 
reference. In this case three different reference frames are involved. Let S denote the local 
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FIG. 1: Illustration of the reference frames involved when comparing the mean free path A g 
measured in frame £ to the formation time of the emitted gluon r" = l/k± measured in frame 
Pi, p- and p'l are the momenta of the incoming particles 1 and 2 in the respective frame, the thick 
dashed arrow (labeled k±) depicts the radiated gluon. See text for more details. 
In this example \p2\ = 2 \pi\ and < (pi,j>2) = 45° are chosen, leading to /3' ~ 0.933 and 6 ~ 69°. 
The gluon in this example is emitted with cosh y = 7" = y/2. 

frame that is co-moving with the average velocity of the medium in each computational cell. 
In this frame the mean free path A 9 is computed from the interaction rates. E' is the center 
of momentum frame of the colliding particles in which the matrix element ([3]) is computed. 
Finally, E" is the reference frame in which the gluon is emitted purely transversal with 
respect to the axis denned by the colliding particles in the CM frame and thus r" = l/k±. 

In order to compare A g to r" in the Theta function modeling the LPM cut-off via 
9(A s -r)=e^- t" \ , the overall boost 

coshy 



7 = T y 1 + p'p» 



[1 + (3' tanhy cos( 



(4) 



from E to E" needs to be taken into account. 7' and 0' denote the boost and the boost 
velocity respectively from E to E'. 7" = coshy and /3" = tanhy are the boost and boost 
velocity from E' to E". The latter can be expressed in terms of the rapidity y of the emitted 
gluon measured from the CM frame E'. 9 is the angle (0 < 9 < 71 / '2) between /?' and the 
axis of the colliding particles in the CM frame as seen from E. See Fig. [TJfor an exemplary 
illustration. 

With this the Theta function entering the bremsstrahlung matrix element can be written 



as 



e ( A 9 -r) = e(^-£) 



cosh y 



tanhy cos( 



(5) 



For thermal energies the boost velocity 0' becomes small,_7 ~ cosh 7/, and the Theta function 
effectively reduces to Q(k±A g — coshy) as employed in [16]. 

For large jet energies the boost factor (j3J reduces the total cross section with respect 
to the simpler expression Q(k^A g — coshy) since the phase space for the k± integration 
gets reduced. More essential for the kinematics of the outgoing particles, however, is the 
peculiar way in which fl5]) distorts the shape of the available phase space. Most notably 
the rapidity of the emitted gluon in the center of momentum frame with respect to the 
incoming jet momentum p[ (CM) gets strongly shifted to negative values with increasing 
7, see lower panel of Fig. [2J For comparison the upper panel of Fig. [2] shows the rapidity 
Ay = arctanh(/3') associated with the boost from laboratory to CM frame. Note that y and 
Ay are not additive due to the angle 9 in 01]). 

While for thermal energies, 7 ~ coshy, the available phase space for the rapidity y is 
essentially given by a 1/ coshy shape and thus symmetric around y — 0, for larger jet 
energies the boost velocity /?' becomes large and the emission in the CM frame is strongly 
shifted to the backward direction. With this, even for small k±, the energy of the emitted 
gluon can become large in the CM frame but will still be small in the laboratory frame due 
to the boost. 

Note that this is still within the scope of the approximations underlying the Gunion- 
Bertsch matrix element [lsj]. Besides the requirement that the transverse momenta are 
small, q±,k± <C y/s, the main approximation that leads to the result ((3]) is that xq± <C k±, 
where x is the fraction of the light cone "+"-momentum of the jet particle carried away by 
the radiated gluon. More specifically, in light cone coordinates, k = (xy/s,qj_/(xy/s),q±,0) 
describes the kinematics of the radiated particle. The limit xq± <C k± includes x ~ k±/y/s 
corresponding to emission into the central rapidity region, y « 0, that is the main interest 
of the authors in But it also includes the x < k±/ y/s region corresponding to backward 
emission, y < 0, and thus the shift towards negative rapidities caused by our effective 
implementation of the LPM effect via a cut-off does not break the validity of the result by 
Gunion and Bertsch that is employed in (J3]). 
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FIG. 2: (Color online) Upper panel: Distribution of the magnitude of the boost from the laboratory 
frame to the center of momentum frame (CM) for different jet energies E (laboratory system) and 
T = 400 MeV expressed in terms of the rapidity Ay = arctanh(/3'). 

Lower panel: Rapidity distribution of the emitted gluon in the CM frame with respect to the 
incoming gluon jet momentum p{ (CM) for different jet energies E (laboratory system). The 
medium temperature is T = 400 MeV. ' 



III. JETS IN A STATIC MEDIUM WITH gg o gg INTERACTIONS ONLY 



To gain a better understanding of the underlying mechanisms of jet quenching in heavy 
ion collisions, it is important to study the evolution of jets in a simplified setup. For 
this we track the evolution of a high-energy gluon as it propagates through a static and 
thermal medium of gluons. To begin with, this section focuses on BAMPS calculations 
exclusively incorporating binary gg -H- gg interactions in order to provide a baseline for 
further investigations. 

The straightforward way to implement such a setup within a microscopic transport model 
such as BAMPS would be to populate a static box with gluons according to a thermal 
distribution at a given temperature T, then to inject a jet particle with initial energy E 
and to track its propagation through the medium. To cut down on computation time, 
however, we choose a more direct, Monte Carlo-type approach. For a jet particle with given 
energy E a certain number of collision partners is generated from a thermal distribution 
with temperature T without actually propagating any medium constituents. This method 
neglects possible effects of the propagating jet on the medium. To ensure consistency, we 
have checked this approach against full calculations of static systems within BAMPS, where 
the dynamics of all particles are explicitly simulated. 

The mean energy loss per unit path length dE/dx is then computed as follows (c = 1) 

where i denotes the interaction type (gg — > gg, gg — > ggg and ggg — > gg) and R 1 is the 
interaction rate for process i. (AE % ) is the mean energy loss in a single collision of type i 
computed as the weighted sum 

The individual weighting factor P- is proportional to the interaction probability for a given 
process. For gg — > gg and g_g^ — > ggg it is Pj = a l jV j , for gluon annihilation processes 



999 ~^ 99 h is P.- = E I 2? E 16]. Here, a 1 is the cross section for a given interaction type, 



E 1 E 2 E 3 

v Tel = 2eIe 2 ^ s ^ e relative velocity of two incoming massless particles and J32 is the phase 
space integral over the matrix element |A4 999 _> g9 | 2 . For the purpose of this section only 
99 ~^ 99 interactions are considered. We will drop this restriction in section IIVI 
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FIG. 3: (Color online) Differential energy loss (dE/dx)22 of a gluon jet in a static and thermal 
medium of gluons with temperatures T = 300, 400 and 500 MeV taking only binary gg — y gg 
interactions into account. 

Figure [3] shows the mean differential energy loss dE/dx of a gluon jet caused by binary 
gg — y gg interactions as a function of the jet energy E and for different medium temperatures. 
The energy loss is computed as described above and the medium is represented by a thermal 
ensemble of gluons. The jet particle is traced such that after each collision the outgoing 
particle with the highest energy is considered the new jet particle. 

The differential energy loss exhibits the expected (see [h]] for a concise overview) log- 
arithmic dependence on the jet energy E and the dominant quadratic dependence on the 
medium temperature T 

dR _ _ / A F!T\ 

(8) 



dE 

dx 



ocC R na 2 s T 2 \n (-§H , 

2^2 V m D J 



where Cr is the quadratic Casimir of the propagating jet, Cr = Ca = N c for gluons. For 
T = 400 MeV and a jet energy E = 50 GeV we find an elastic energy loss of ^\ 2 _^ 2 = 
1.2 GeV/fm that increases to %\ 2 ^ 2 = 2 GeV/fm at E = 400 GeV. 

More detailed information than in the mean energy loss per unit path length is contained 
in the time evolution of the energy distribution of the jet particle propagating through the 
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medium. As for the differential energy loss a Monte Carlo approach is chosen where the 
collision partners are sampled from a thermal distribution. Given discretized time steps 
At, the gluon density n(T) and a fixed number of collision partners N, the implemented 
interactions are sampled according to their probabilities 

n(t)At n(t)At I 32 n(t) 2 At 

P 22 = v rel a 22 ^- Fa-v^-jr- Pja = ^-^-^—. (9) 

Starting at t = Ofm/c with an initial parton energy Eq this approach yields the evolution of 
the energy as a function of time, E(t). Repeating this procedure many times we compute 
p(E) dE, the probability that a parton that started with E(t = Ofm/c) = Eq has an energy 
E < E(t) < E + dE at a given time t. 

Figure H] shows the probability distribution for the energy of a gluonic jet particle injected 
with an initial energy of Eq = 50 GeV after certain propagation times in a thermal medium 
of T = 400 MeV and T = 600 MeV respectively. In both cases the distribution of the jet 
energy induced by collisions with the constituents of the medium becomes rather broad. A 
distinct peak at lower energies only re-emerges at very late times, roughly after 50fm/c for 
T = 400 MeV and 30fm/c for T = 600 MeV. The mean energy loss as depicted in Fig. [3] 
is therefore a valuable observable but contains only limited information. It is noteworthy 
that there exists a finite probability for the jet to gain energy by collisions with the thermal 
gluons. This effect is more pronounced for higher medium temperatures. 

As already found in I4J] the shapes of the distributions induced by collisional energy loss 
significantly differ from models that employ a mean energy loss accompanied by momentum 



diffusion such as 
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IV. JETS IN A STATIC MEDIUM INCLUDING INELASTIC gg o ggg INTER- 
ACTIONS 

Having investigated the evolution of a gluon jet that exclusively undergoes elastic gg — > gg 
interactions, we now make use of the full potential offered by BAMPS and include inelastic 
99 ^ 999 processes according to the matrix element ([3]). As in the previous section we 
explore the propagation of gluonic jets in a static and thermal medium of gluons to system- 
atically analyze the energy loss mechanism in simulations with the microscopic transport 
model BAMPS. As for the purely elastic case, the observables presented in this section are 
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FIG. 4: (Color online) Time evolution of the energy distribution of a gluon jet that interacts only 
via gg — > gg processes with a static and thermal medium of gluons with T = 400 MeV (upper panel) 
and T = 600 MeV (lower panel). The initial (t = Ofm/c) energy of the gluon jet is Eq = 50 GeV. 
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computed by means of the Monte Carlo approach presented in section II III 
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FIG. 5: (Color online) Differential energy loss of a gluon jet in a static and thermal medium of 
gluons with T = 400 MeV. The contributions of the different pQCD processes implemented in 
BAMPS {gg —> gg, gg —> ggg and ggg —> gg) to the total dE/dx are shown. 



Figure [5] shows the mean differential energy loss dE/dx of a gluon jet with energy E in a 
static thermal medium with T = 400 MeV caused by binary gg — > gg and inelastic gg O ggg 
interactions. The contributions from the different processes to the total energy loss are 
displayed separately. From this compilation it is obvious that bremsstrahlung processes gg — > 
ggg are by far the most dominant contribution to the gluonic energy loss within the BAMPS 
framework, whereas gluon annihilation processes are negligible and binary interactions as 
discussed in detail in section II I II contribute only on a small level. The resulting differential 
energy loss is almost linearly rising with the energy, for example resulting in a total dE/dx ~ 
32.6 GeV/fm at E = 50 GeV. 

At large energies the temperature dependence of the resulting total differential energy 
loss appears to be linear as can be seen from Fig. [6j where dE/dx of a gluon jet is compared 
for T = 400 MeV, T = 500 MeV and T = 600 MeV. This stems from the dominant gg — > ggg 
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FIG. 6: (Color online) Total differential energy loss, including gg f-> ggg processes, of a gluon jet 
in a static and thermal medium of gluons at T = 400 MeV, T = 500 MeV and T = 600 MeV. 

processes and is in contrast to the elastic energy loss ([H]) that exhibits a quadratic dependence 
on the temperature. Possible logarithmic contributions to the temperature dependence 
cannot be resolved at this stage. 

Given these results it is necessary to discuss the origin of the strong energy loss from 
radiative processes within BAMPS. First of all, despite the large differential energy loss for 
gluonic jets, the individual cross sections increase only slowly with the jet energy as seen 
in Fig. [7J For instance the average total cross sections for a gluon jet with E = 50 GeV 
in a gluonic thermal bath with a temperature T = 400 MeV are (a gg -^gg) ~ 1.3 mb and 
(<Jg g ^ggg) ~ 3.5 mb. This emphasizes that BAMPS does indeed operate with reasonable 
partonic cross sections based on pQCD matrix elements. 

Since the cross section for gg — > ggg processes as seen in Fig. [7J yields moderate mean 
free paths for jet-like particles, the cause for the large differential energy loss needs to be a 
large mean energy loss per single collision, (AE). The energy carried away by the radiated 
gluons, however, is in itself not sufficient to explain such large mean AE. From the radiation 
spectrum of a gluon jet with E = 50 GeV, shown in the upper panel of Fig [8] for different 
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FIG. 7: (Color online) Cross section of a gluon jet with energy E interacting with thermal gluons 
(T = 400 MeV). Shown are the cross sections for gg — > gg processes and for gg — > ggg processes. 

medium temperatures, a mean energy of the radiated gluon can be read off that is below 
the (AE) that would be needed to fully explain the magnitude of dE/dx, a finding we 
will explicitly confirm later. The spectrum is displayed as the number of radiated gluons 
per energy interval du and per distance dx scaled by the total number of emitted gluons 
and weighted with the gluon energy. Here u is the lab frame energy of the gluon that 
in the center of momentum frame is emitted with transverse momentum k± according to 
the Gunion-Bertsch matrix element ([3]). The spectra are clearly peaked at energies that are 
small compared to the energy of the parent jet, with a small tail reaching out to high energies. 
With increasing temperature the peak of the spectrum shifts towards higher energies in an 
apparently linear way, favoring the emission of gluons with higher energies. 

For completeness, the lower panel of Fig. [8] shows the angular distribution of gluons 
radiated off a E = 50 GeV gluon jet in a T = 400 GeV medium for different ranges in the 
energy of the radiated gluon u. The angle ( is taken in the lab frame with respect to the 
initial direction of the parent jet. With increasing energy u the radiated gluons are emitted 
more preferably at small angles, only for soft gluons there is a sizable probability to be 
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FIG. 8: (Color online) Upper panel: Energy spectrum ^ N ^ dx of radiated gluons per energy 
interval dco and distance dx. w is the energy (lab frame) associated with the radiated gluon 
according to the Gunion-Bertsch matrix element ([3j). The radiating jet is a gluon with E = 50 GeV 
traversing a thermal medium with T = 300 MeV, T = 400 MeV or T = 500 MeV. 
Lower panel: Angular distribution of the radiated gluon in the lab frame with respect to the 
original jet direction for different energies u of the radiated gluon. Jet energy E = 50 GeV, 
medium T = 400 MeV. 
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emitted transversely or in the backward direction. However, as is clearly visible in Fig. [BJ 
due to the cut-off in transverse momentum (j3J) the gluons cannot be emitted at very forward 
angles, an effect that is more pronounced for low u>. 

To emphasize that the energy carried away by the radiated gluon alone does not explain 
the observed mean energy loss per collisions, it is noteworthy that the strong and linear rise 
in the energy loss due to gg — > ggg is only present when identifying the outgoing particle with 
the highest energy as the outgoing jet and thus using AE = E m — max (E° nt , E% nt , E^ nt ). 
This is the most natural choice and is employed throughout all calculations in this work. 
The average energy oj of the radiated gluon, however, is rising much slower with the jet 
energy, as can be seen in the upper panel of Fig. |HJ This is due to the fact that the energy 
is distributed among three outgoing particles, the gluon emitted with energy u being only 
one of them. In fact, assuming u < E m /3 (in agreement with the results presented in Fig. 
[9]), AE min = co is only the smallest possible energy loss, while the largest energy loss allowed 
by energy and momentum conservation is AE max = E m — ( ) = Ein + U . 



Looking at the distribution of AE in the lower panel of Fig. [9] for a fixed jet energy 
that underlies the mean energy loss per collision as shown in the upper panel of Fig. [9j it 
becomes obvious that it is indeed a very fat tail of the AE distribution that causes a large 
averaged (AE). The distinct peak at low AE can be readily identified as being related to 
energy carried away by the radiated gluon. 

The reason for the heavy tail in the AE distribution needs to be looked for in the compli- 
cated plethora of configurations for the outgoing particles allowed by the underlying matrix 
element ([3]). In the appendix, figures IT9l and 1201 show some examples for such configurations 
that were randomly chosen according to the matrix element ([3]). The most specific feature 
of the phase space sampled in gg — > ggg is that the radiated gluon is predominantly emitted 
into the backward hemisphere in the center of momentum frame, compare © and Fig. [2j 
Due to the strong bias towards negative rapidities that is present for large boosts (jl]), the 
energy of the radiated gluons in the CM frame is in many cases comparable to the energies 
of the two other outgoing particles even for small transverse momenta k±. 

With this some typical scenarios can be identified. Fig. [18] in the appendix |X] illustrates 
that for a given cut on the momentum transfer qj_ there are distinct configurations correlating 
the energy of the radiated gluon in the center of momentum frame with the energy of one of 
the remaining outgoing particles. Note that due to the nature of the Gunion-Bertsch matrix 
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FIG. 9: (Color online) Upper panel: Energy loss AE23 of a gluon jet in a single gg —> ggg process. 
Two different cases are shown: One where AE = E m — max (E° ut , E% nt , E% ut ) (the upper group 
of lines) and one (labeled "rad.", the lower group of lines) where AE = uj is the energy of the 
radiated gluon. 

Lower panel: Distribution of AE 23 , with AE = E in - max (E% nt , E% ut , E$ nt ) , for E = 400 GeV 
and T = 0.4 GeV. 17 



element ([3]) the transverse momentum k± is typically comparable to the momentum transfer 
q±. For example, selecting on configurations where q± and E[ are small and E' 3 , i.e. the 
energy of the radiated gluon in the CM frame, is large, yields an energy loss significantly 
above its most probable value, thus these configurations contribute to the heavy tail observed 
in Fig. HI This is especially true for cases where particle 1 is emitted into the forward 
hemisphere, because for these specific configurations the outgoing particle 2 needs to be 
emitted into the forward hemisphere as well in order to guarantee energy and momentum 
conservation. Boosted back into the lab frame the available energy is thus mainly split 
between the particles 1 and 2 yielding a large energy loss. The same line of reasoning holds 
for cases where E[ is on the order of E' 3 , see example Fig. I19dl in the appendix IBl 

This illustrates that the fat tail in the AE distribution is mainly caused by configurations 
where in the center of momentum frame the radiated gluon is emitted with a large energy into 
the backward hemisphere and the remaining energy is split among the two other particles 
going into the forward hemisphere. Events with large q± and k± (see Fig. Il9bl and Fig. 
I20al in the appendix for examples) also yield a large energy loss but are significantly less 
probable due to the steeply falling l/q\ contribution in the matrix element (J3]). Compare 
appendix [A] and especially Table [I] for a more quantitative analysis of the above mentioned 
configurations. 

As already discussed in connection with the scenario containing only elastic interactions, 
valuable information beyond the mean energy loss is contained in the evolution of the energy 
distribution p(E) dE of the jet. Starting out with p(E) = 5 (E — Eq) we have seen in 
Fig. H] that elastic collisions cause a broadening of the distribution with a distinct peak 
at low energies only re-emerging at very large times. Because of the much stronger mean 
energy loss caused by inelastic gg — > ggg processes a more rapid evolution is to be expected 
when all interactions included in BAMPS are taken into account. Indeed, Fig. [10] shows 
that the energy distribution of a jet with Eq = 50 GeV traversing a gluonic medium with 
T = 400 MeV is spread over almost the entire range after roughly lfm/c. A distinct peak 
at£si 7.5T emerges at about 3.5 fm/c for T = 400 MeV and 2 fm/c for T = 600 MeV. 

A quantity that is often used to characterize the effect of the medium on a jet-like particle 
is q. It is defined as the sum of the transverse momentum transfers squared divided by the 
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FIG. 10: (Color online) Time evolution of the energy distribution of a gluon jet in a static and 
thermal medium of gluons with T = 400 MeV (upper panel) and T = 600 MeV (lower panel). The 
initial (t = Ofm/c) energy of the gluon jet is Eq = 50 GeV. Inelastic gg -H- ggg interactions are 
taken into account. 
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path length L the particle has traveled 



(10) 



where i runs over all collisions the particle has undergone within the path length L. From 
now on we will simply use the time t instead of L since we are dealing with massless particles. 
Alternatively, if one knows the average momentum transfer squared per mean free path as 
a function of the jet energy one can compute the mean q as 



Typically q is used to quantify the transverse momentum picked up from elastic collisions 
that eventually induce the radiation of Bremsstrahlung gluons. In the commonly used 
eikonal approximation the jet particle acquires no additional transverse momentum due to 
the radiation of the Bremsstrahlung gluons. In our approach however, radiative and elastic 
interactions are treated on equal grounds and jets can also pick up transverse momentum in 
inelastic gg — > ggg processes. In the following, we therefore naturally extend the definition 
of q as given above to also describe the evolution of transverse momentum due to inelastic 
gluon multiplication processes within BAMPS. 

As a cross-check we have compared the result from both approaches, ( flOl) and (TlTl) . 
using independent calculations and found perfect agreement. The upper panel of Fig. [TT1 
shows the average momentum transfer squared per mean free path (Ap 2 ± )/X as a function 
of the jet energy in a gluonic medium with T = 400 MeV. A logarithmic behavior at 
large energies can be seen for (Ap"j_)/X from binary gg — > gg interactions, with {Ap\)/X « 
2.3 GeV 2 /fm at E = 50 GeV rising to {Ap\)/X « 3.7 GeV 2 /fm at E = 400 GeV. As reflected 
in the differential energy loss, the average transverse momentum transfer squared per mean 
free path for inelastic gg — > ggg interactions is much higher, (Ap±)/X ~ 22.8GeV 2 /fm at 
E = 50 GeV and (Ap\)/X rj 64.2GeV 2 /fm at E = 400 GeV while the gluon annihilation 
processes ggg — > gg virtually do not contribute at all. 

The lower panel of Fig. [TT] shows (q) as defined in equations (TTU|) and (ITT]) as a function 
of the path length L = t for a gluon jet with initial energy E = 50 GeV. As before the 
medium is characterized by T = 400 MeV. Over the range up to t = 3.5fm/c shown in Fig. 
[TTJ the contribution from elastic interactions is almost constant at (5)22 ~ 2.3GeV 2 /fm. 
For jets that interact only via binary gg — > gg one finds that (g) 2 2 is actually slowly and 




(11) 
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linearly falling to (9)22 ~ 1.9GeV 2 /fm at t — 50fm/c. The combined (q) is dominated 
by the radiative gg — > ggg contribution and starts at (q) w 23GeV 2 /fm, falling to (q) ~ 
12.5GeV 2 /fm at t = 3.5fm/c. This indicates that the negligence of transverse momentum 
pick-up in radiative processes might indeed be an oversimplification. 

The numbers for q found in this work are well within the range of values found by 
other theoretical energy loss schemes, though the comparison is difficult since q in these 
calculations often is a free parameter or related to free parameters. Fitting to experimental 



data the authors of 



13] have found go f° r the central region of Au+Au at r = 0.6fm/c, 



where conditions should be roughly comparable to the setup used in this section, to be 
ranging from 2.3 GeV 2 /fm based on the Higher Twist approach 22[, over 4.1 GeV 2 /fm based 

to 18.5GeV 2 /fm based on the 



up 



24] . In 25] the application of the Higher 



on the approach by Arnold-Moore- Yaffe (AMY) 
approach by Armesto-Salgado- Wiedemann (ASW) / 
Twist approach to jet quenching data yields go ~ 3.2GeV 2 /fm for a medium evolution 
based on BAMPS (employing r = 0.3fm/c), while a hydro based medium evolution yields 
g ~ 0.9GeV 2 /fm (r = 0.6fm/c). 



V. NON CENTRAL AU+AU COLLISIONS AT RHIC ENERGY 



In a previous work [15] BAMPS has been applied to simulate elliptic flow and jet quench- 
ing at RHIC energies, for the first time using a consistent and fully pQCD-based microscopic 
parton transport model to approach both ke y o bservables of heavy ion collisions at RHIC 



within a common setup. As established in [26|, [27] the medium simulated in the parton 
cascade BAMPS exhibits a sizable degree of elliptic flow in agreement with experimental 



findings at RHIC, while showing a small ratio of shear viscosity to entropy r]/s [28]. The 
suppression of high-p^ gluon jets in central, 6 = fm, collisions has been found to be roughly 
constant at Raa^ ~ 0.053 with radiative events gg — > ggg being the dominant contribution 
to jet energy loss. This nuclear modification factor is in reasonable agreement with recent 



analytic results for the gluonic contribution to the nuclear Raa [1Q(, though the suppression 
of gluon jets in BAMPS appears to be slightly stronger. However, we expect improved agree- 
ment in future studies when employing a carefully averaged (b) that will be better suited 
for comparisons to experimental data, having a centrality range even for the most central 
selections, than the strict b = fm case. 
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FIG. 11: (Color online) Upper panel: Average transverse momentum transfer squared per mean 
free path, (Ap^_)/A, for gg — > gg, gg — > ggg, ggg — > gg and for the sum of all processes. The 
medium temperature is T = 400 MeV. 

Lower panel: (q) (t) as defined in (fit)]) and (fTTI) for gg — > gg processes and for all BAMPS processes 

as a function of the path length L = t for a gluon jet with initial energy Eq = 50 GeV. The jet 

99 

traverses a medium as characterized above. 



In order to further test the energy loss as modeled by BAMPS and to allow for more 
extensive comparisons with experimental data and analytic models, our investigations need 
to be extended. One way to gain further insight is to study the evolution of jets in the 
simulations of Au+Au collisions including light quark degrees of freedom and to subsequently 
employ a fragmentation scheme that would allow for direct comparison with hadronic data. 
The consistent implementation of inelastic 2- and 3-particle interactions for all possible 
combinations of light quarks, light anti-quarks and gluons is challenging and we leave this 
to an upcoming work. Another way is to look at more differential observables, for instance to 
study the nuclear modification factor as a function of centrality or to study the dependence 



on the angle with respect to the reaction plane [29|, |30 |. 

Here we will present first results in this direction by discussing the nuclear modification 
factor as simulated by BAMPS for non-central Au + Au collisions at the RHIC energy of 
y/s = 200AGeV with a fixed impact parameter b = 7fm. This roughly corresponds to 20% 
to 30% experimental centrality. 

For these simulations the initial gluon distributions are sampled according to a mini-jet 
model with a lower momentum cut-off po = 1.4 GeV and a i^-factor of 2. The underlying 
nucleon-nucleon collisions follow a Glauber-model with a Wood-Saxon density profile and 
the results by Gliick, Reya and Vogt [32 1 are used as parton distribution functions. Quarks 
are discarded after sampling the initial parton distribution since currently a purely gluonic 
medium is considered. To model the freeze out of the simulated fireball, free streaming 
is applied to regions where the local energy density has dropped below a critical energy 
density e c . This setup has been successfully checked against experimental findings such 
as the distribution of transverse energy in rapidity and the flow parameter u 2 at various 



centralities in 



17,126, 2 



The test particle method [16| is employed to ensure sufficient statistics and to allow for 
the resolution of adequate spatial length scales. For the calculations of the b = 7fm events 
presented here, we use iVtest = 220. The size of the computational grid in the transverse 
direction is Ax = Ay = 0.2 fm and in the longitudinal direction an adaptive grid linked to 
the number of particles in the cells is used that roughly corresponds to Arj ~ 0.1 -v- 0.2. 

Because of the steeply falling momentum spectrum of the initial gluon distribution the 
computational expense of calculating observables at high-p^ becomes very high, rendering a 
brute force approach that consists in simulating a huge number of random events completely 
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FIG. 12: (Color online) Gluonic nuclear modification factor Raa at mid-rapidity (ye [—0.5, 0.5]) 
as extracted from simulations of Au+Au collisions at 200 AGeV with fixed impact parameters 
b = Ofm and b = 7fm. For comparison experimental results from PHENIX 
for central (0% - 10%) and off-central (20% - 30%) collisions. 

infeasible. Therefore, in an approach similar to the concept of importance sampling, only se- 
lected events are chosen for computation and a suitable weighting and reconstruction scheme 
is applied. A large number of initial spectra is characterized according to X = max ((pr)J, 
the maximum p? of the particles in an event in a given rapidity range. Subsequently X is 
divided into bins of appropriate size and from each bin j in X a number Nj of events is 
selected for simulation. The results are averaged over the Nj events within a bin and finally 
combined with the appropriate weights Pj, where Pj is the probability that in a simulated 
initial state there are one or more particles with their transverse momentum in the interval 
(pr) 3 -_i > ipr)j but none with pt > (pt)j- We use a bin size of Apr = 1 GeV and for this 
work select up to Nj = 40 events per bin. 

We have successfully tested this method by reconstructing simple analytic distribu- 
tions such as p(s) = e~ s and more importantly by reconstructing the initial mini-jet pt~ 
distribution from a small set of selected events. 
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Since within the mini-jet model employed in these simulations the scaled nucleon-nucleon 
Pr _ spectrum is directly accessible, the usual definition of the nuclear modification factor 

R AA = d2N AA/dp T dy 
AA T AA d 2 a NN /dp T dy 

simplifies to just the ratio of the final p^-spectrum over the initial p^-spectrum. 

Figure [12] shows the gluonic contribution to the nuclear modification factor R AA as ex- 
tracted from simulations with BAMPS for fixed impact parameters b = Ofm and b = 7fm. 
The critical freeze out energy density is e c = 1.0GeV/fm 3 , see section IVT1 for results em- 
ploying e c = 0.6GeV/fm 3 . Most notably the quenching pattern for b = 7fm is flat as for 
the central collisions but at R AA ~ 0.13. This is qualitatively consistent with experimental 
observations since the measured shape of the nuclear modification factor for 7r° does not 
change significantly going from the centrality class 0% - 10% to 20% - 30% [311 ] . 

A comparison in terms of the magnitude of the jet suppression for b = 7fm is difficult 
since to the best of our knowledge there are no published results from analytic models that 
explicitly disentangle the gluon and quark contributions to non-central R AA . Taking the 
ratio of non-central to central R AA presented in Fig. [12] as a rough guess indicates that the 
decrease in quenching is more pronounced in BAMPS than in the experimental data. The 
ratio of the nuclear modification factor between central (0% - 10%) and more peripheral 
(20% - 30%) collisions is Raa\o%-io% I -^^l2o%-30% ~ 0.6 f° r ^ ne experimental data, while 
for the BAMPS results RAA\b=of m / -R-AA\b=7f m ~ 0-4- However, for thorough quantitative 
comparison one needs to address the question how fragmentation influences the relative 
gluon and quark contribution to the nuclear modification factor at different centralities. 
The issue of detailed quantitative comparison therefore needs to be settled once light quarks 
and a fragmentation scheme are included into the simulations. 

To complement the investigations of Raa at a non-zero impact parameter, in Fig. [13] 
we present the elliptic flow parameter t> 2 for gluons computed within BAMPS at the same 
impact parameter b = 7fm as used above. The elliptic flow in this figure is obtained from 
simulations using a critical energy density e c = 0.6GeV/fm 3 in order to be comparable 
to previous results. Going from e = lGeV/fm 3 to e — 0.6GeV/fm 3 does not change the 
magnitude of the nuclear suppression factor R AA but slightly distorts its flatness, see the 
results marked "X = 1" in Fig. [TT] Compared to previous studies the elliptic 

flow results presented here extend to high-p^ gluons, up to transverse momenta of roughly 
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FIG. 13: (Color online) Elliptic flow parameter v<i for gluons as a function of transverse momentum 
pr- Extracted from BAMPS calculations of a non-central b = 7fm Au + Au collisions at 200 AGeV 
using e c = 0.6GeV/fm 3 . Shown are two cases: The regular scenario where X = 1 and the scenario 
in which the LPM cut-off is modified by a factor X = 2. 



Pt ~ 10 GeV. Bearing in mind the limited statistics of this calculation, the Vi of high- 
er gluons is rising up to pr ~ 4 GeV. Afterward, from about pt ~ 5 GeV on, the results 
indicate a decrease in the elliptic flow for high-p^ gluons. This behavior is in good qualitative 
agreement with recent RHIC data [33] that for charged hadrons shows v 2 to be rising up to 
t> 2 ~ 0.15 at px ~ 3 GeV followed by a slight decrease. 



VI. SENSITIVITY OF THE RESULTS ON THE LPM CUT OFF 

Every effective model comes with a set of free parameters whose choice can be motivated 
by physical arguments or by fits to experimental data. The most notable free parameter 
in the transport model BAMPS is the coupling strength a s that has been fixed to the 
canonical value of a s = 0.3 throughout this work. The consequences of differei 
a s , especially on elliptic flow observables, have been studied in previous works 
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simulations of heavy ion collisions the freeze out energy density e c , see section [V] is also a 
parameter that can be adjusted within certain limits. 



4 



X = 0.5 
X= 1.0 




1.5 - 



1 



0.5 - 







50 



100 



150 



200 



250 



300 



350 



400 



E [GeV] 



FIG. 14: (Color online) Total cross section for gg — > ggg processes involving a gluon jet with energy 
E that traverses a thermal medium (T = 400 MeV). Shown are three different values for the factor 
X that modifies the effective LPM cut-off (fl3|) . 

When investigating inelastic gg <H- ggg interactions within BAMPS there basically enters 
another parameter due to the effective modeling of the LPM effect via a cut-off. In the 
99 ~ >* 999 matrix element ([3]) the cut-off is realized via a Theta function that essentially 
compares the formation time r of the radiated gluon to the mean free path of the the parent 
gluon A g as discussed in section [III in more detail. However, the argument underlying the 
distinction between coherent and incoherent processes via a threshold r = A g is of course a 
qualitative one. When effectively modeling the LPM effect via a cut-off the Theta function 
fl5]) could therefore be replaced by a more general form 




(13) 



where X is a real number not too far from 1. 
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In this section we explore the consequences on our results when modifying the LPM- 
cutoff ()5]) by a factor X. Specifically we choose X = 0.5, X = 1 (the usual choice) and 
X = 2. This should provide a grasp on how sensitive the results for partonic energy loss 
and collective flow within BAMPS simulations are on the specific prescription for including 
the LPM effect. 

The computation of total cross section for radiative gg — > ggg processes involves an inte- 
gral of the matrix element ([3]) over the transverse momentum k± . It is therefore straightfor- 
ward that a larger X in the cut-off (fl3l) corresponds to a smaller total cross section. Though 
the actual dependence is not linear and indeed non-trivial, Fig. HH that shows the total 
cross section for gg — > ggg processes in a gluonic medium with T = 400 MeV for different 
choices of X, confirms the simple qualitative considerations. The change in the cross section 
naturally corresponds to a change in the rate for this process via R23 = (^o"23) and thus in 
the mean free path between radiative processes A 2 3 = 1/-R23- 



250 




E [GeV] 

FIG. 15: (Color online) Differential energy loss of a gluon jet with energy E that traverses a 
thermal medium (T = 400 MeV). Shown are three different values for the factor X that modifies 
the effective LPM cut-off ([13]) . 
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FIG. 16: (Color online) Energy loss AE23 of a gluon jet in a single gg — > ggg process (T = 400 MeV) . 
Two different cases are shown. One where AE = E m — max [E° nt , £v? ut , -Kj Ut ) (the upper group 
of lines) and one (labeled "rad. only", the lower group of lines) where AE = uj is the energy of 
the radiated gluon. For each case three values for the parameter X in the LPM cut-off (|13p are 
explored. 

Correspondingly the differential energy loss dE/dx is affected by a change in the param- 
eter X as shown in Fig. [T^J Larger X leads to a larger mean free path and thus a smaller 
energy loss per path length. That the change in dE/dx is indeed mainly due to the change 
in the total cross section can be seen when comparing with Fig. [161 where it is shown that 
the effect of changes in X on the energy lost in a single gg — > ggg interaction is rather small. 
Also note that the ordering of the energy loss in a single interaction is not as straightforward 
as in the total cross section. This is due to the highly non-trivial impact of the LPM cut-off 
on the sampling of outgoing particle momenta in gg — > ggg processes. 

Finally we investigate the effect of a modified LPM cut-off f lT3|) on observables in full 
simulations of heavy ion collisions as discussed in section [V] Due to limited computing 
resources we restrict ourselves to a comparison of the cases X = 1 and X = 2. As was to 
be expected from the change in the energy loss in a static medium (Fig. [13]) comparing 
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FIG. 17: (Color online) Comparison of gluonic Raa for different values of the parameter X in the 
effective LPM cut-off (|13p . Raa extracted at mid-rapidity (y e [—0.5, 0.5]) for Au+Au collisions 
at 200 AGeV with fixed impact parameter b = 7fm. Critical freeze out energy density is e c = 
0.6GeV/fm 3 . 

different cut-off parameters, the level of jet quenching is considerably reduced when going 
from X = 1 to X = 2, see Fig. [TTJ Note that in the calculations for Fig. [T7] we have used 
a freeze out energy density of e c = 0.6GeV/fm 3 , which leads to a slight linear decrease of 
Raa towards high p? as compared to the results with e c = 1.0GeV/fm 3 in Fig. [T2j The 
elliptic flow as seen in Fig. [T3]is reduced by roughly 30% to 35% when going to X = 2, but 
the qualitative features as a function of pt remain unaffected. 

VII. SUMMARY 

In this work we have presented detailed results on the energy loss mechanisms for high- 
energy gluons within the microscopic parton transport model BAMPS. To this end we 
have investigated the evolution of high energy gluons within thermal and static media of 
gluons. This setup is known as the "brick problem" and is proposed by the TECHQM 
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collaboration as means of comparing results from different parton cascade models and Monte 
Carlo calculations. We have discussed collisional energy loss as well as the contribution of 
radiative processes implemented in BAMPS via the Gunion-Bertsch matrix element Q. The 
inelastic gg — > ggg processes are found to be the dominant source of energy loss for high 
energy gluons in computations within the BAMPS model resulting in a strong differential 
energy loss that rises almost linearly with the jet energy. The strong mean energy loss in 
99 ~~ 999 processes is due to a heavy tail in the AE distribution for single interactions, 
caused by the phase space configurations of outgoing particles dictated by the Gunion- 
Bertsch matrix element ([3]) in combination with the effective LPM cutoff This prefers 
the emission of radiated gluons into the backward hemisphere with energies that in the center 
of momentum frame are comparable to that of the remaining outgoing particles, while they 
are small in the laboratory frame. The jet energy in these cases can be split into two particles 
yielding a large energy loss. 

For a purely gluonic medium with T = 400 MeV the transport parameter q as defined in 
(flOl) and (fill stemming from binary gg — > gg interactions is found to be roughly constant 
at q = 2.3GeV 2 /fm. When including inelastic gg -H- ggg processes, q as a measure of the 
accumulated transverse momentum exhibits a stronger dependence on the path length and 
is much larger than for elastic interactions, q = 12 -j- 23 GeV 2 /fm. 

Extending our investigations of jet quenching in central Au + Au collisions at RHIC 
energies |l5j |. we have presented the gluonic contribution to the nuclear modification factor 
for non-central b = 7fm collisions simulated in BAMPS. The gluonic Raa in these off- 
central events is found to be flat at Raa ~ 0.13 over a wide range in for a critical energy 
density of e c = 1.0GeV/fm 3 , being in qualitative agreement with the experimental results. 
Results employing e c = 0.6GeV/fm 3 exhibit a slight linear decrease in Raa towards high 
transverse momenta. 

Since our transport approach allows for the simultaneous investigation of high-pr observ- 
ables and bulk properties of the medium, we have also studied the elliptic flow parameter 
v 2 for gluons up to roughly 10 GeV. v% peaks at a pr ~ 4 5 GeV and slowly drops towards 
larger transverse momenta. 

We have explored the sensitivity of the results on our implementation of the LPM effect 
by varying the momentum space cut-off that enters all computations involving the matrix el- 
ement Mg g ^.g gg for the inelastic processes included in BAMPS. Changing this cut-off mainly 
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affects the total cross section for radiative gg — > ggg processes, leading to changes in the 
energy loss per mean free path. The nuclear modification factor is found to be more sensitive 
on the specific implementation of this cut-of than the elliptic flow. 

The characteristics of the strongly interacting, but still fully pQCD based, medium within 
the BAMPS description will be studied in further detail and in improved quantitative com- 
parison to experimental data in upcoming works including light quark degrees of freedom. 
This of course needs to be accompanied by a hadronization scheme that allows for a more 
direct connection to hadronic observables. While a consistent modeling of low-p^ hadroniza- 
tion needs careful consideration, the application of fragmentation functions to the high-p^ 
sector will be straightforward once light quarks are included. 

The application of the transport model BAMPS to heavy quark, i.e. charm and bottom, 
elliptic flow and quenching will provide further valuable insight into the underlying mech- 
anisms and the importance of a careful treatment of the fireball dynamics. Such studies 
are underway and will complement the results presented in this work. Recently an imple- 
mentation of radiative processes based on Gunion-Bertsch type matrix elements has been 
successfully applied to heavy quark t>2 and Raa 34], indicating promising perspectives for 
the implementation within the established BAMPS framework. 

Furthermore, in future projects we will study the medium response to high-p^ particles 
that deposit energy into the medium created in heavy ion collisions via energy loss mecha- 
nisms. This will for example allow us to investigate the possible formation of so called mach 
cones. First studies 35] have already very successfully demonstrated that BAMPS offers 
the ability to describe collective shock phenomena in a viscous hydrodynamic medium. 
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Appendix A: Typical phase space configurations in gg —> ggg processes 

The characteristics of the kinematics in gg — > ggg processes dictated by the Gunion- 
Bertsch matrix element in combination with the LPM cut-off can be studied by identi- 
fying typical regions of the phase space and by selecting on these regions. Respecting energy 
and momentum conservation, the kinematics of the three outgoing particles can be described 
by 6 independent parameters. One such possible choice would be the set (E[, E 3 , cos(#i), 
cos(6 , 3 ), 0i, 3 ), where E[, E' 3 are the energies, 6\ and #3 are the angles with respect to the 
incoming momentum pj et and <pi and <fr 3 are the azimuthal angles of the outgoing particles 1 
and 3 (the emitted gluon), with all values being taken in the CM frame. The corresponding 
values for the outgoing particle 2 can then be inferred from momentum conservation. 

Another choice would be to replace cos(#i) and cos(#3) by the momentum transfers q± 
and k±_ as directly given in fl3]). Note however, that this choice hides the information whether 
cos(#i) and cos(#3) are larger or smaller than zero, i.e. whether particles 1 and 3 are emitted 
in the forward or in the backward direction. Finally, also replacing E 3 by y, the rapidity of 
the emitted gluon, would yield a set of parameters that are closest to matrix element (E]). 

For the purpose of this discussion we use the set (E[, E' 3 , q^, k±, l5 <p 3 ) and throw in 
additional information by looking at the signs of cos(#i) and cos(#3) as needed. We ignore 
any dependence on the azimuthal angles <pi, <p 3 and note that for a fixed q± the matrix 
element fl3]) gives a k± that is typically on the order of q±. Thus we select events according 
to E[, E 3 and q±. 

Considering an E = 400 GeV jet-like gluon inside a thermal medium with T = 400 MeV, 
we select events having a low q±, GeV < q± < 3 GeV, and events having a rather high q±, 
8 GeV < q± < 12 GeV. Figure [TBI then shows the color coded correlations between E[ and 

For low q±, GeV < q± < 3 GeV, two distinct regions in the E^-.E'3-plane are visible. 
The energy of the emitted gluon is quite high in all cases due to the strong preference of 
events with y < caused by the LPM cut-off as discussed in section [TTJ One region features 
small E[, comparable to q± while the other regions features large E[ f» E 3 . For large q±, 
8 GeV < q± < 12 GeV only one distinct region emerges with both E[ ~ E' 3 large. 

Table UJ lists the mean energy loss for events within these kinematical regions, differen- 
tiating between emission of particle 1 into the forward direction (cos(#i) > 0) and into the 
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< q± < 3 


< E[ < 5 


12 < E' 3 < 20 


cos(0i) < 
cos(6>i) > 
-1 < cos(6»i) < 1 


7.5% 
7.5% 
15.0% 


{AE) ps 13GeV 
(AE) ps 61 GeV 
(AE 1 ) ps 37 GeV 


12 < E[ < 20 


12 < E 3 < 20 


cos(0i) < 
cos(6»i) > 
-1 < cos(6»i) < 1 


0.4% 
13.6 % 
14.0 % 


(AE 1 ) ps 7.7 GeV 
(AE 1 ) ps 30.7 GeV 
{AE} ps 30.0 GeV 




8 < q± < 12 


8 < E( < 15 


8 < E' 3 < 15 


cos(6»i) < 
cos(6>i) > 
-1 < cos(6»i) < 1 


0.34 % 
0.56 % 
0.9% 


(AE) ps 105 GeV 
(AE) ps 135 GeV 
(AE) ps 123 GeV 



TABLE I: Mean energy loss for given cuts in q±, E[ and E' 3 , cf. Fig. [181 The percentage given in 
column 4 corresponds to the fraction of all events within these cuts relative to the total number of 
events. For reasons of readability the unit statement GeV is omitted for q±, E[ and Eg. 

backward direction (cos(#i) < 0). Additionally the relative abundance of events within these 
regions relative to all events is given. 
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FIG. 18: (Color online) Probability distribution of the outgoing energies E[ and Eg in the center 
of momentum frame (CM) for gg —> ggg processes given a certain cut in the momentum transfer 
qj_, where E' 3 is the energy of the emitted gluon. Jet energy E = 400 GeV, medium temperature 
T = 400 Me V. 

Upper panel: GeV < q± < 3 GeV. Lower panel: 8 GeV < q± < 12 GeV. 
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Appendix B: Illustration of randomly selected gg — > ggg events 

In order to visualize the possible configurations of outgoing particles in gg — > ggg events, 
the figures [TH] and [2U] illustrate a number of gg — > ggg events in the center of momentum 
frame. These events feature an incoming E = 400 GeV gluon jet that interacts with con- 
stituents from a thermal medium (T = 400 MeV). These events have been randomly selected 
obeying the relative importance given by the matrix element ([3]). All events are rotated such 
that the incoming jet momentum in the center of momentum frame points along the positive 
x-direction and that the outgoing momentum p[ is in the x-y-plane. 



36 




(a) AE = 0.67 GeV, 

cos(Gi) = 0.998, p{ = (22.3, 22.2, 1.4, 0.0) 
cos(8 2 ) = -0.487, p' 2 = (0.2, -0.1,0.2,0.0) 
cos(9 3 ) = -0.997, p' z = (22.2, -22.1, -1.6, -0.0) 
(E[ ah , E l 2 ah , E\ ah ) = (399.3, 1.0, 1.4) GeV 
q±_ = 1.43 GeV, k±_ = 1.60 GeV 



(b) AE = 206.43 GeV, 

cos(6i) = -0.134, p[ = (30.2,-4.0,29.9,-0.0) 
cos(9 2 ) = 0.829, p' 2 = (5.3, 4.4, 2.2, 2.0) 
cos(6 3 ) = -0.011, = (32.1,-0.3,-32.1,-2.0) 
{E[ ah ,E l 2 ah , £ 3 ab ) = (152.5, 56.9, 193.6) GeV 
qx = 29.89 GeV, k x = 32.13 GeV 



Pz 




(c) AE = 26.01 GeV, (d) AE = 117.01 GeV, 

cos(Gi) = -0.365, p{ = (4.3, -1.6, 4.0, 0.0) cos(Gi) = 0.998, p{ = (23.0, 23.0, 1.5, -0.0) 

cos(e 2 ) = 1.000, p' 2 = (22.8, 22.8, 0.3, 0.5) cos(6 2 ) = 0.976, p' 2 = (9.4, 9.1, 1.7, 1.1) 

cos(8 3 ) = -0.980, p^ = (21.6, -21.2, -4.3, -0.5) cos(6 3 ) = -0.994, p' z = (32.3, -32.1, -3.3, -1.1) 

(Ef° , E\ ah , Ef°) = (21.8, 374.0, 6.0) GeV {E[ ah , E l 2 ah , £ 3 ab ) = (283.0, 112.7, 9.2) GeV 

q±_ = 4.02 GeV, k ± = 4.34 GeV q x = 1.54 GeV, kj_ = 3.45 GeV 

FIG. 19: (Color online) Randomly selected gg — > ggg events involving a gluon jet with E = 
400 GeV (lab system) displayed in the CM frame (T = 400 MeV). All events are rotated such 
that the incoming jet momentum (CM) points along the positive x-direction and that the outgoing 
momentum p[ is in the x-y-plane. All kinemsSrcal values are given in the CM frame, execpt for 
the energies of the outgoing particles in the lab frame, I?J ab , E^, -Ei ab = ui. Part 1: Events 1-4 




(a) AE = 54.26 GeV, 

cos(Gi) = 0.861, p{ = (31.5, 27.1, 16.0, -0.0) 
cos(e 2 ) = 0.983, p' 2 = (3.0, 3.0, 0.5, 0.2) 
cos(9 3 ) = -0.876, vk = (34.4, -30.1, -16.5, -0.2) 
{E[ ah , E l 2 ah , E\ ah ) = (345.7, 35.2, 22.6) GeV 
q± = 16.03 GeV, k±_ = 16.54 GeV 



(b) AE = 71.64 GeV, 

cos(6i) = 0.997, p{ = {17 A, 17 A, 1.2, -0.0) 
cos(6 2 ) = -0.277, p' 2 = (8.0, -2.2,-1.2, 7.6) 
cos(8 3 ) = -0.894, = (17.0, -15.2, -0.0, -7.6) 
{E[ ah ,E 2 ah ,El ab ) = (328.4, 50.9, 22.2) GeV 
q±_ = 1.25 GeV, k±_ = 7.60 GeV 




(c) AE = 14.15 GeV, 

cos(Gi) = -0.133, p{ = (1.4, -0.2, 1.3, 0.0) 
cos(e 2 ) = 0.999, p' 2 = (16.6, 16.6, -0.1, 0.7) 
cos(8 3 ) = -0.996, p^ = (16.5, -16.4,-1.2, -0.7) 
(E[ ah , E l 2 ah , E l f h ) = (14.0, 385.8, 1.0) GeV 
q±_ = 1.35 GeV, k ± = 1.39 GeV 

FIG. 20: (Color online) Same 



(d) AE = 41.35 GeV, 
cos(Gi) = 0.953, p{ = (24.1, 23.0, 7.3, 0.0) 
cos(e 2 ) = -0.881, = (22.3,-19.6,-9.3,4.9) 
cos(6 3 ) = -0.535, p^ = (6.3,-3.4,2.0,-4.9) 
{E[ ah ,E 2 ah ,El ab ) = (358.7, 17.3, 26.8) GeV 
q x = 7.28 GeV, k±_ = 5.34 GeV 

.[H Part 2: Events 5-8 out of 8. 
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